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ABSTRACT 

Coalescence of super massive black holes (SMBH’s) in galaxy mergers is potentially the d ominant 
contr ib utor to the low frequency gravitational wave background (GWB). It was proposed by Merritt & 


Ekers (2002) that X-shaped radio galaxies are signposts of such coalescences, and that their ab un dance 
might be used to predict the magnitude of the gravitational wave background. In Roberts et al. (20151 
we present r adio images of all 52 X-shaped radio source candidates out of the sample of 100 selected 
by Cheung (2007) for which archival VLA data were available. These images indicate that at most 
21% of the candidates might be genuine X-shaped radio sources that were formed by a restarting of 
beams in a new direction following a major merger. This suggests that fewer than 1.3% of extended 
radio sources appear to be cand idates for genuine axis reorientations (“spin flips”), much smaller than 
the 7% suggested by|Leahy & Parma ((1992). Thus the associated gravitational wave background may 
be substantially smaller than previous estimates. These results can be used to normalize detailed 
calculations of the SMBH coalescence rate and the GWB. 

Subject headings: galaxies: active — gravitational waves — radio continuum: galaxies 


1. INTRODUCTION 

The gravitational wave background from SMBH coa¬ 
lescence during galaxy mergers is proving to be a diag¬ 
nostic of galaxy formation models because of the useful 
constraints on the GWB from im proved und ersta nding 
of timing noise in pulsar data (Shannon et al. J2013). Re¬ 
cent studies of the expected rate of merger-based SMBH 
formation based on current observations (corrections to 
BH-bulge relation, and new est im ates of g alaxy merger 
timescales) (|McWilliams, Ostriker, fe Pretorius 2013| 


Sesana |2013a| |Ravi et al.||2014a|) suggest that the rates 
derived earlier from semi-analytic models for gal axy for¬ 
mation implemented in dark matter simulations (Sesana, 


Vecchio, & Colacino 
be seriously underes 


2008] |Sesana & Vecchio||2010|) may 
imatedC This has led to the predic¬ 


tion of imminent detection of the GWB by pulsar tim¬ 


ing arrays (Siemens et al. 2013). On the other hand, 
the GWB from coalescing supermassive binaries might 


tions 

Sesana 

2013b 

Ravi et al. 

2014b| and owing to 

the “last parsec" pro! 

)lem (Milosavljevic & Merritt 2003 

Colpi 2014 

Vasiliev| 2014 

). 


tain inversion symmetric morphologies in radio galaxies 
where presence of close companion galaxy and binary 
SMBH models ha ve been used in some cases to explain 


the phenomenon (Ekers et al. 1978 Be gelman, Bland- 
ford, fe Rees||1980| |Hunstead et al~ |1984| |Taylor et aT* 

19901 IDunn, Fabian, fc Sanders|2006 lErlund et al. 120061 


Subrah manya n~et al. 2006| Saripalli et al. 2013). Mer¬ 

ritt & Ekers ( 2002| suggested that the structures of radio 
galaxies, which are hosted by giant ellipticals, offer inde¬ 
pendent and direct clues to the prevalence of SMBH bina¬ 
ries. In particular, they pointed out that “X-shaped” ra¬ 
dio galaxies (XRG’s) may be signposts of SMBH coales- 

roberts@brandeis.edu 


cence and hence their abundance amongst the extended 
radio source population might be an indicator of the true 
coalescence rate. Using the estimate of the f raction of 
radio galaxies that show X morphology from |Leahy fc] 


Parma 


(1992) and a relic lobe lifetime of 10 8 yr they es- 
timated the merger rate to be about 1 Gyr -1 galaxy -1 . 

X-shaped radio sources in particular, and the direc¬ 
tional stability of the jets in active galactic nuclei in gen- 


binary formation (e.g., Ekers et al. (1978 

): Begelman, 

Blandford, & Rees 

(1980); Gower et al. ( 

1980); Deane 


context depends on their formation mechanism - whether 
their radio structures arise from reorientations of the BH 
spin axes (spin flips), or slow precession or drifts in the 
axes, or whether the X-shaped structures are caused by 
the hydrodynamic deflection of backflows or jet-shell in¬ 
teractions or have nothing to do with changing direc¬ 
tions of the central beams (for references an d a discus- 
sion of the various proposed mechanisms see [Roberts et| 


al. (20151. These different mechanisms may be operating 
in different cases, generating the variety in radio struc¬ 
tures; therefore, it is crucial to determine the specific 
abundance of “genuine XRG’s” that are candidates for 
reorientations of SMBH spin axes as these are what are 
relevant for the GWB. All of the above mechanisms for 
generating off-axis distortions lead to radio structures 
with low axial ratios, and to deduce the precise mech¬ 
anism in individual cases might be difficult. For this 
reason, good quality radio, optical, and X-ray studies 
of populations of complete samples of extended radio 
sources that have small axial ratios are critical to the 
problem. 


In Roberts et al. (20151 we presented an initial step 
in this direction using available data in the Very Large 
Array archi ves on the sam ple of 100 XRG candidates 
compiled by Cheung ( 2007) from the NRAO FIRST sur- 
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vey. Characterization of the extended structure in all of 
the sources showed that a large majority of the sample 33 
of 52 sources, or 63%) constitutes sources where the off- 
axis emission is traced to individual lobes. These sources 
were further divided into those where the emission occurs 
near the inner ends of the lobes (25 of 52 sources, or 48%) 
and those where the emission occurs near the outer ends 
(8 of 52 sources, or 15%). 

An example of each of these classes is shown in Fig¬ 
ure [TJ The sources with distortions near the inner ends 
of the on-axis lobes, and in which the off-axis emission 
on the two sides are not collincar, appear to be the best 
candidates for back flow deflections by therma l halos of 
the host ellipticals ( Leah y & Williams (1984); but see 
Gong, Li, & Zhang (201lJ). It is possible that such inver¬ 


sion symmetric structures (which would classify as XRGs 
in the literature if the off-axis emission is long enough) 
might arise from rapid swings or flips of the jet axis dur¬ 
ing a merger; however, in explaining the variety in source 
structures this would require models where in some cases 
nuclei along with the host galaxy are displaced a con¬ 
siderable distance during mergers. Additionally, in most 
cases the timescales associated with the flip would require 
to be fine tuned so that on one hand it is short enough 
to produce distinct wings and lobes and long enough for 
the older lobes of edge-brightened type to have morphed 
into edge-darkened or FR-I type. Hydrodynamic forces 
associated with the thermal gas halo of the host ellip¬ 
tical would now have to be invoked to create the emis¬ 
sion gap between the two lobe pairs. Most importantly, 
wings are almost always along the host minor axis and 


the active main lo bes a re towards the major axis (Capetti 
et al. 2002 Saripalli & Subrahmanyan 2009), also for a 


sub-sa mple of the present set of sources, |Roberts et al.| 
(2015)): this requires that all of these inversion symmet- 
ric XRGs would have to be formed as a result of mergers 
that flip the axis from the minor to major axis. 

On the other hand, the sources with distortions near 
the outer ends of the on-axis lobes appear to be candi¬ 
dates for cases of “axis precession or drift,” with radio 
structur es t h at may signal the presence of binary black 


holes (Begclman, Blandford, & Rees 1980 Deane et al. 


20141, and thus be precursors of genuine X systems. Fx- 


cluding these two source categories, eleven sources (21%) 
have been identified as genuine X-shaped radio galaxy 
candidates Q Two examples of candidate “genuine X” 
sources are shown in Figure [2] 

We point out that our source morphology classifica¬ 
tion into outer-end and inner-end deviation sources can 
have implications for the general population of sources 
presently classified as XRGs. In our classification scheme 
several of the well-known XRGs will be classified as inner- 
end deviation sources and as argued above are more eas¬ 
ily explained via the backflow deflection phenomenon 
rather than the spin-flip model. For the sources in Fig¬ 
ure [2] the observations that we have do not indicate that 
the central transverse emission originates in either of the 
two lobes. Hence they are classified as X-shaped radio 
sources. 

In the present paper we discuss the implications of 


1 These sources are: J0144-0830, J1008+0030, J1015+5944, 
J1043+3131, J1327—0203, J1345+5233, J1406+0657, J1408+0225, 
J1606+0000, .11614+2817 and J1625+2705. 


50 02 47 

i i i i i i i i i 

f 17 


c 

o 


Declinati 

o 

ro 

iilL ■' 

01 47 




32 

- (a) 


EL_ i _i_i_i_i_i_i_i_ 


07 02 52 51 50 49 48 47 46 45 44 

Right ascension (J2000) 



08 45 11 10 09 08 07 06 

Right ascension (J2000) 


Fig. 1.— VLA L-band images of the radio galaxies 
(a) J0702+5002 (lowest contour = 0.2 mJy/beam, peak = 
7.17 mJy/beam), and (b) J0845+4031 (lowest contour = 
0.15 mJy/beam, peak = 20.8 mJy/beam). These are examples 
of sources where the off-axis emission occurs (a) at the inner ends 
of the main lobes and (b) at the outer ends of the lobes. Note the 
prominent inversion symmetry in each source’s structure. 

these results for the GWB. 


2. GRAVITATIONAL RADIATION FROM SMBH 
COALESCENCE 

2.1. Electromagnetic Signatures of Axis Reorientations 

Our imaging of a sample of low-axial ratio radio sources 
has allowed us to identify examples where the central 
transverse emission may be genuinely independent ra¬ 
dio emission unrelated to the current activity. These 
are sources in which, within the limitations of sensitiv¬ 
ity and resolution of the observations, the central trans¬ 
verse emission on the two sides cannot be traced to either 
of the two lobes, and appears as a swathe of emission 
through the central core or galaxy but at a large angle 
with respect to the main radio axis. We have identified 
11 such sources; these sources are out of the sample of 52 
sources whose radio structures we have examined care¬ 
fully. These 52 are out of t he sample o f 100 sources that 
have been identified by Cheung (2007) to have low axial 
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Fig. 2.— VLA L band images of radio galaxies that are can¬ 
didates for genuine X-shape sources created by axis reorienta¬ 
tions. (a) The upper panel shows J1043+3131 (lowest contour 
= 0.3 mjy/beam, peak = 44.6 mjy/beam); (b) the lower panel 
shows J1406+0657 (lowest contour = 0.15 mjy/beam, peak = 
71.4 mjy/beam). 


ratios and hence represent sources with substantial off 
axis emission. The sample itself was selected from the 
FIRST survey that was searched for extended sources by 
selecting those with fitted major axis exceeding 15" (the 
FIRST beam was about 5" FWHM); this yielded a total 
of 1648 extended radio sources. The search of the FIRST 
images was limited to sources with peak flux density ex¬ 
ceeding 5 mJy beam -1 . Altogether this leads to an esti¬ 
mate that at most 1.3% of extended radio sources might 
be genuine X-shaped radio sources. If in these sources 
(which we call genuine XRG’s) the central swathe of off- 
axis emission arises because of axis reorientation alone 
then our observations point to an occurrence rate of less 
than 1.3% for axis reorientations amongst the extended 
radio source population. 

It may be noted that this rate of axis-reorientations 
in radio galaxies is much smaller than the occurr ence 
rate (of 7%) estimated by Leahy & Parma (1992) for 
XRGs (in the luminosity range of 0.3-30 x 10" W Hz -1 ). 


However, Leahy & Parma derived only the abundance 
fraction of radio galaxies with wings, without examin¬ 
ing whether they were genuine X-shaped radio sources 
and without excluding sources where the off-axis emis¬ 
sions were obviously inversion-symmetric extensions of 
the lobes of a double lobed structure. Based on the 
FIRST survey, the number of s uch cand idate X-shaped 
sources was deemed to be 100 (Cheung 2007), which is 
indeed a 6% fraction of the 1648 extended radio sources 
that were selected for examination. 

On a related note, if the eight sources that appear to 
have off-axis emission connected to the outer lobe ends 
are representative of axis precession or drifts then our 
observations point to an occurrence rate of 0.9% for axis 
drifts (occurring over a smaller angle of few tens of de¬ 
grees). 

It has been argued that sources with substantial off- 
axis emission that is traceable to the individual inner- 



sources, which also included these sources with wings 
traceable to the inner-lobe ends as well, have a statisti¬ 
cal relationship with the shape of the host elliptical in 
that the current radio axis is aligned closer to the major 
axis of the galaxy. We see no obvious reason why axis 
reorientations ought to be almost always from an initial 
direction along the host minor axis to one along the ma¬ 
jor axis. It is more reasonable therefore to suppose that 
the off-axis emission in most sources in such an inclusive 
sample arises from ba ckflow d efle ctions, as concluded by 


Saripalli & Subrahmanyan (2009), and that the number 
of genuine X-shaped sources is a small subset of these 
candidates, as we suggest here. 


2.2. Gravitational Wave Background from SMBH 
Coalescence 

Prediction of the level of the GWB arising from the 
mergers of m assive galaxies depend s on estimation of the 


merger rate (Jaffe & Backer 2003). This has been done 
in a vari ety of ways including studies of close pairs of 


Xu et al. 2012 Jiang, Jing, & Han 2014| 

Keenan et al. 

2014|, simulations of structure formation ( 

Kitz bidder & 

White 2008; Guo et al. 2011: Rodriguez-"Gomez et al. 


ber of close galaxy pairs ( Keenan et al.|2014 ) show that a 


typical L* galaxy has undergone ~ 0.2 —0.8 mergers since 
z ~ 1. Taking the mean of these numbers implies a ma¬ 
jor merger rate somewhat less than 0.1 Gyr -1 galaxy -1 . 
Estimation of the level of the GWB also depends on un¬ 
derstanding the process of SM BH co alescence as d riven 
by dynamical friction |^hmgyj£ Tremaine 19871 and 


how the last parsec problem 

is overcome 

Milosavljevic 

& Merritt 

2003 

Colpi]2014 

Vasiliev 2014 

. It would be 


tion for all of this, and the X-shaped radio galaxies offer 
one avenue in that direction. 


2.3. Constraints Arising from the XRG Population 
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In this model, the coalescence rate for luminous ellip¬ 
ticals might be estimated by examining the occurrence 
rate of genuine X-shaped radio sources in the population 
of extended radio sources. If we use our finding above 
that at most 1.3% of extended sources might have gen¬ 
uine X-shapes, then the occurrence rate of coalescence 
in their hosts is given by 0.013/r, where r is the lifetime 
of relic lobes of radio galaxies. If we adopt the value of 
r = 1 0 8 yr for th e lifetimes of the relics (as assumed by 
Merritt & Ekers (2002)) then we arrive at an estimate 
of 0.13 Gyr _i galaxy -1 as an upper limit to the event 
rate for radio galaxies, and presumably therefore for the 
parent population of massive ellipticalsr] On the other 
hand, if we include the eight sources with off-axis emis¬ 
sion traced to outer ends of lobes as sources that have 
undergon e precessi on due to the presence of binary black 
holes ( Deane et al.p 014) then we estimate a coalescence 
rate of at most twice this value, of ~ 0.2 Gyr -1 galaxy -1 . 
In either case the event rate is substantially smaller than 
the event rate of about 1 Gyr -1 galaxy -1 usually assumed 
in calculations of t he GWB from galaxy formation (Mer¬ 
ritt fe Ekers||2002 1. 

Our revision of the event rate is owing to the mor¬ 
phology classification adopted. If we consider the four 
examples of XRGs in [Merrit t & Ekers (2002) in their 
Figure 2, only 3C 223.1 would have been included as a 
potential “genuine XRG” as per our scheme. For the 
other three XRGs (3C52, 3C403, & NGC326) given the 
available images we would note them as cases where the 
wings link with individual lobes rather than forming an 
independent transverse emission feature. 


3. SUMMARY AND CONCLUSIONS 


In Roberts et al. (20151 we analyzed 1.4 and 5 GHz 
archival VLA continuum data on a sample of 52 FIRST 
radio sources selected on the basis of low-axial ratio ra¬ 
dio structures. In the current paper we discuss the im¬ 
pactions of the results for the gravitational wave back¬ 
ground. Our principal results are: 


1. We argue that a large fraction of the sample (63%) 
consists of sources where the off-axis emission is 
traced to individual lobes, which are more likely to 
be due to back flows or axis drifts rather than axis 
flips. 


2. In this case the number of genuine X-shaped radio 
sources is at most eleven (21%). 


3. We then derive a five times lower occurrence rate 
for genuine X-shaped radio ga l axies than the value 
estimated by Leahy & Parma (1992). 


4. If X-shaped radio galaxies arise because of major 
mergers alone then our observations point to an 
occurrence rate of less than 1.3% for such mergers 
amongst the extended radio source population. 


5. The low occurrence rate for major mergers sug¬ 
gests a significantly smaller event rate for coales¬ 
cences of supermassive black holes of less than 


2 The radiation generated by SMBH coalescence in t he few 
larges t galaxies is expected to dominate the GWB (Sesana et al. 
2004}. 


0.13 Gyr -1 galaxy -1 for hosts of radio galaxies. 
Inclusion of radio galaxies with off-axis emission 
traced to outer ends of lobes as sources that have 
undergone axis precession due to a close binary 
black hole only increases the coalescence rate by 
a factor of two. 

6. If the class of sources with inner deviations are also 
cases of axis flips, then the event rate would be 
substantially greater, and similar to that derived by 
Merritt and Ekers (2002). Indeed a measurement 
of the GWB might be a clue to the origins of the 
XRG structures. 

7. Our estimate of the coalescence rate of giant ellip¬ 
ticals could be used to normalize calculations of the 
GWB from such events. 


4. FUTURE INVESTIGATIONS 

The elucidation of the structures in non-classical radio 
sources is important in that it has the power to yield 
event rates for SMBH coalescence and occurrence rates 
for a variety of phenomena that perturb the stability of 
the BH spin axis and interrupt the accretion flow that 
power s the beams. The work presented in |Roberts et al.| 
(20151 and herein using primarily VLA archival data has 


been impactful in its conclusions for the event rate that 
determines the GWB, and clearly the conclusions could 
be substantially firmer with better data. Improved imag¬ 
ing in radio and optical broad band continuum, followed 
by derivation of source orientations in 3D with respect 
to the host ellipticity using depolarization distributions, 
and deriving the abundance statistics for the different 
types of off-axis distortions, may reveal the formation 
mechanisms for small-axial-ratio radio sources and hence 
provide a useful constraint on GWB from BH coalescence 
that accompanies interactions that reorient the spins of 
radio loud AGNs. 

While the improved resolution (compared to th e 
FIRST surv ey) in the observations presented in Roberts] 
et al. (|2015) has helped in unraveling basic structure m 
several”!’ the sources, the sensitivity is too limited in 
tracing fully the larger lobe and wing structures, all of 
which are critical in isolating genuine X-shaped radio 
galaxies. Deeper observations are needed to determine 
the structures more accurately than what have been im¬ 
aged with the limited sensitivity archival observations 
presented here. Radio observations at multiple wave¬ 
lengths would be useful in detecting gradients in spectral 
index within the lobes and hence indicate the time evolu¬ 
tion in the structure. We anticipate that with improved 
data the number of genuine candidates for axis reorienta¬ 
tions would be further constrained, thereby reducing the 
event rates further and the discordance with formation 
simulations would be on firmer ground. 

Critical to the estimation of the SMBH binary forma¬ 
tion rate is the assumed lifetime of the relic radiation in 
an XRG. The estimate of 10 8 yr ( Merritt fe Ekers||2002 1 
was based primarily on modeling of the spectral aging 
of the o ff -axis wings. With improved observations of the 
Cheung (2007) sample we will be able to make similar 
estimates for a, larger sample of XRGs and thus firm up 
this critical parameter in the computation of the coa¬ 
lescence rate. Thus new JVLA observations are needed 
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to test different formation models, estimate occurrence 
rate of axis reorientations, and hence also to place more 
accurate limits on the expected GWB signal. 

If models for axis flips are the explanations for inner 
lobe distortions, then we may explore the consequence 
of the inference that flips almost always occur from mi¬ 
nor to major axis. Indeed, giant radio sour ces do have 


their axes alonj 

the 

Subrahmanyan 

2009) 


plorations of the torques and tidal couplings between the 
large-scale structure of the host elliptical and the inner 


accretions disk and hence BH spin axis may suggest an 
explanation for this phenomenology. 


5. ACKNOWLEDGMENTS 

The National Radio Astronomy Observatory is a fa¬ 
cility of the National Science Foundation, operated un¬ 
der cooperative agreement by Associated Universities, 
Inc. D. H. R. gratefully acknowledges the support of the 
William R. Kenan, Jr. Charitable Trust. We thank the 
referee, Prof. David Merritt, for his useful suggestions. 


REFERENCES 


Begelman, M. C., Blandford, R. D., Sz Rees, M. J. 1980, Nature, 
287, 307 

Binney, J., & Tremaine, S. 1987, Galactic Dynamics, (Princeton, 
NJ: Princeton Univ. Press), 747 
Capetti, A., Zamfir, S., Rossi, P., Bodo, G., Zanni, C., Sz 
Massaglia, S. 2002, A&A, 394, 39 
Carlberg, R. G., Cohen, J. G., Patton, D. R., et al. 2000, ApJ, 
532, LI 

Cheung, C. C. 2007, AJ, 133, 2097 
Colpi, M. 2014, http://arxiv.org/abs/1407.3102 
Deane, R. P., Paragi, Z., Sz Jarvis, M. J. et al. 2014, Nature, 511, 
57 

Dunn, R. J. H., Fabian, A. C., Sz Sanders, J. S. 2006, MNRAS, 
366, 758 

Ekers, R. D., Fanti, R. Lari, C., Sz Parma, P. 1978, Nature, 276, 
588. 

Erlund, M. C., Fabian,, A. C., Blundell, K. M., Celotti, A., Sz 
Crawford, C. C. 2006, MNRAS, 371, 29 
Gong, G. P., Li, Y. P. and Zhang, H. C. 2011, Ap.J, 734, L32 
Gower, A. C., Gregory, P. C., Unruh, W. G., Sz Hutchings, J. B. 
1982, ApJ, 262, 478 

Guo, Q., White, S., Boylan-Kolshin, M. et al. 2011, MNRAS, 413, 
101 

Hunstead, R. W., Murdoch, H. S., Condon, J. J., Sz Phillips, 

M. M. 1984, MNRAS, 207, 55 
Jaffee, A. H., Sz Backer, D. C. 2003, ApJ, 583, 616 
Jiang, C. Y., Jing, Y. P., & Han, J. 2014, ApJ, 790, 7 
Keenan, R. C., Focaud, S., DePropris, R. et al. 2014, ApJ, 795, 
157 

Kitzbichler, M. G. & White, S. D. M. 2008, MNRAS, 391,1489 
Leahy, J. P., Sz Williams, A. G. 1984, MNRAS, 210, 929 
Leahy, J. P., Sz Parma, P. 1992, in Extragalactic Radio 

Sources. From Beams to Jets, ed. J. Roland et al. (Cambridge: 
Cambridge Univ. Press), 307 
McWilliams, S. T., Ostriker, J. P., Sz Pretorius, F. 2013, 
http://arxiv.org/abs/1211.5377 
Merritt, D. Sz Ekers, R. D., 2002, Science, 297, 1310 


Milosavljevic, M., Sz Merritt, D. 2003, in The Astrophysics of 
Gravitational Wave Sources, AIP Conf. Proc., 686, 201 
Ravi, V., Wyithe, J. S. B., Shannon, R. M., Sz Hobbs, G. 2015, 
MNRAS, 447, 2772 

Ravi, V., Wyithe, J. S. B., Shannon, R. M., Hobbs, G., Sz 
Manchester, R. N. 2014, MNRAS, 442, 56 
Rodriguez-Gomez, V., Genel, S., Vogelsberger, M., 2015, 

MNRAS, 449, 49 

Roberts, D. H., Cohen, J., P., Liu, J., Saripalli, L., Sz 
Subrahmanyan, R. 2015, ApJS, in press 
Saripalli, L., Malarecki, J. M., Subrahmanyan, R., Jones, D. H., 

Sz Staveley-Smith, L. 2013, MNRAS, 436, 690 
Saripalli, L., Sz Subrahmanyan, R. 2009, ApJ, 695, 156 
Sesana, A. 2013a, MNRAS, 433, LI 
Sesana, A. 2013b, CQGra, 30, 4014 

Sesana, A., Haardt, F., Madau, P., Volonteri, M. 2004, ApJ, 611, 
623 

Sesana, A., Sz Vecchio, A. 2010, Class. Quant. Grav., 27, 084016 
Sesana, A., Vecchio, A., Sz Colacino, C. N. 2008, MNRAS, 390, 
192 

Shannon, R. M., Ravi, V., Coles, W. A., et al. 2013, Science, 342, 
334 

Siemens, X., Ellis, J., Jenet, F., Romano, J.D. 2013, CQGra, 30, 
4015 

Subrahmanyan, R., Hunstead, R. W., Cox, N. L. J., Sz McIntyre, 
V. 2006, ApJ, 636, 172 

Taylor, G. B., Perley, R. A., Inoue, M., Kato, T., Tabara, H., Sz 
Aizu, K. 1990, ApJ, 360, 41. 

Thorne, K. S. 1987, in Three Hundred Years of Gravitation, ed. 

S. W. Hawking Sz W. Israel (Cambridge: Cambridge Univ. 
Press), 330 

Vasiliev, E. 2014, http://arxiv/ab/1411.1762vl 
Wen, Z. L., Liu, F. S., Sz Han, J. L. 2009, ApJ, 692, 511 
Xu, C. K., Zhao, Y., Scoville, N., Capak, P., Drory, N., Sz Gao, Y. 
2012, ApJ, 747, 85 




